REPLICA SYMMETRY BREAKING IN AN AXIAL MODEL OF 

QUADRUPOLAR GLASS 

N. V. Gribova, E.E.Tareyeva 

We perform the replica symmetry breaking (RSB) in the vicinity of the point 
of instabihty of the rephca symmetric solution in the model of axial quadrupolar 
glass. It is shown that the solution with the first stage RSB is stable against the 
second stage RSB. Although there is no reflection symmetry the IRSB solution 
bifurcates continuously from the RS one. These facts mean that our model can not 
be associated with one of two classes usually considered in spin glass theory. 
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1. In recent years the interest to non-Ising spin glasses without reflection sym- 
metry is renewed. The theory of these models is often associated with the theory of 
real structural glasses. At present time there is no satisfactory microscopic model 
of the liquid-glass transition, although a great number of real and computer exper- 
imental data exist, as well as some phemenological theories. At the same time the 
investigations of the spin glasses, orientational glasses and other models with frozen 
disorder are actively carried out. They generated lots of ideas and methods that 
are used for investigation of other systems, for example neural networks, proteins, 
rubbers and real glasses. 

Usually two main links between the spin glass theory and the theory of real 
glasses are marked. First, the spin glass theory of some class of models is used as a 
possible scenarium for the liquid-glass transition in many-particle systems (see, for 
example, Second, the methods of the spin glass theory are used to construct 

a model for glass-forming systems of real particles 

In fact, there is a third link: the approaches are developed where the transition 
to the multipole glass state appears as a part of the liquid-glass transition[|^ |[. 
It is worth to notice that the physical meaning of the order parameters is different 
in different approaches. The problem considered in this paper can be useful in 
connection with the first and the third aspects of links. 

As the recent works have shown, experimental characteristics of relaxation pro- 
cesses are described well enough by the mode coupling theory 0. The similar 
equations can be obtained also as a result of the dynamic investigation of some spin 
glasses. This similarity is most pronounced (the first observation of this type is due 
to the authors of |T^) in the class of mean- field spin glasses without reflection sym- 



metry. In this case the static transition at the temperature Tc (the replica symmetry 
breaking - RSB) is discontinuous and the one-step RSB (IRSB) solution occurs to 
be stable, while the full Parisi scheme fails. The absence of the reflection symmetry 
means that the internal symmetry group does not contain the element transform- 
ing dynamical variable Q to — Q. This results in cubic terms in Ginzburg-Landau 
functional for non-random case ||ll|] and in the special form of RSB free energy func- 
tional for random interactions |T^. The temperature of dynamical transition in 



these models is higher than T^. As a result of the dynamical transition the system 
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is trapped in the state that is less energetically favorable than the RSB state and 
stays there for a long time. The examples of such models are the p-spin model, 
orientational glasses, random Potts model and etc. They can be considered to be 
a prototype of real glasses. Some of these models, especially the p-spin spherical 
model, have already been investigated in detail in the middle of the 90ties (see the 
review [|I|). The others are under investigation today |T^-|T^. It is important 
that the set of the specific features either appears as a whole or is absent as in 
Sherrington-Kirkpatrick model with Parisi solution. The set contains the absence 
of the reflection symmetry, the order parameter discontinuity, the stability of IRSB 
solution and the fact that > Tf.. So traditionally the idea of two classes of models 
is used i, I, p . 



However, there are some indications that this classification is not quite correct. 
For example, the situation is not clear for the three state Potts model ( [p!^ -[p!^). 
The IRSB solution has no discontinuity also in the case of three spin spherical model 
for the certain values of the external field ||2^ . 

In this paper we consider the quadrupolar glass model with infinite-range random 
interaction ETl, 



without reflection symmetry. We show that in contrast to the 
usual classification the IRSB solution is continuous but stable. May be, this is due 
to the influence of the quadrupole operator algebra, that causes internal fields in the 
system with the effect similar to the external field effect in (For more details 
about internal fields in quadrupolar glasses see p3| and [p4| .) 

2. So, we investigate the system of N axial quadrupoles that are in the sites i,j 
of the regular lattice with the Hamiltonian |2! 
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H E JijQiQji 



where Q = 3 — 2, J = 1, = 1, 0, —1, and the values of the coupling constants 
Jij are distributed following the Gaussian law: 

P(J,,) = (v^J)-iexp [-{J,, - J^f/2J% 

and Jo = Jo/N, J = J/N^/\ 

Using the standard procedure of the replica method, we get the expression for 
the free energy, corresponding to the Hamiltonian (|l]) ||2l|, |22| : 



{F)j fj/^^ 1_,„^,., 1 



— lim —max 



NkT n^o n 
J 



+ In Tr exp 



(a/3) 



2 

1/2 



kT ^[kT^ 2{kT)\ 



(2) 



(a/3) 
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Here (a/3) means the sum over the couples of rephcas, n - the number of rephcas, 
and x" ~ {{Q°')t)j, 0°"^ ~ {{Q"Q^)t)j- Further, we suppose for simphcity that 
^0 = 0. 

The case of symmetric rephcas was investigated in details in 



21, 221 and the 



solutions of the equations for the order parameters were obtained. It was shown that 
they grow continuously with the decrease of the temperature. No phase transition 



(in the traditional sense) was found. This fact agrees with the experiment |^5|] for 
N2 — Ar mixture. There was also obtained the behaviour of the specific heat with 
the wide maximum and linear dependence on the temperature for low temperatures, 
that corresponded to the experiment of that time, which had not reached ultralow 
temperatures. 

Later this model was investigated in detail in many papers (see, for example [p^ - 



30[] and references therein) , however the replica symmetry breaking was not carried 
out. However the investigation of the stability of the replica symmetric solution 
using Almeida-Thouless method |^ showed the instability against RSB [^, 

In this work we carry out the replica symmetry breaking for the model (||) in 
the vicinity of the instability point of the RS solution. We find the solution in this 
vicinity, that corresponds to IRSB and we show, that it is stable for further replica 
symmetry breaking. It is important that despite of absence of reflection symmetry, 
IRSB solution bifurcates continuously that contradicts the usual classification. 

To carry out IRSB we divide n replicas into n/m groups with m replicas in each 
group. We take equal to qi if a and (3 belong to one group and q"^ equal to go 
in the opposite case. 



Now: 



E (9'"^) ~ q\n{m — 1) + qon{n — m). 



(3) 



,a/3 



+ 



2ngi + giEQ°- 



(4) 



Substituting (|^) and @ into @, using the expression for linearization of the 
exponent: 



exp(Aa^) 

and using new variables qi 
for the free energy: 



NkT 



1 



dxexp 



-{-mp 
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— + \Pl\ax 



t{p + v), go tp, X 
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kT-> 



(5) 



we obtain 



1) + 4(p + v] 
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- — fdy^ In f dz^m"" {9) 



(6) 



where 



da 



'l-K 



^ = 2e^ + e 



-26 



= tyy/p + tz^ + -{p + v -2 + x), 

and x,p,v,m satisfy the equations that express the extremum condition of the func- 
tional (|6f) : 



X 



dy 



p -\- V 



p 



dy 



J dz^m^^ ' 
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From the expression for free energy (P) we get the entopy in form: 



(10) 



s 

Wk 



-e - -t'x' 
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-mp^ + {p + vY{m- 1) -t'^{p + v) + 2t^x+ 



+ ^Jdy^ In J dz^^!"^ {9). 



3. Expanding (|7|-[10D in t at t 0, it is easy to show, that there is only one 
solution for the system of new equations at high temperatures and it is the same as 
for RS equations [21, 22]. One can obtain the RS equations from (@-10) if f = 0: 



dy^ 

p' = 4: I dy^ 



2e^' + e-2«' 



2e^' + e- 



29' 



:ir 



(12) 



r- V 

9' = ty^p' + -{p' + x' -2). 

The solutions x = x',p = p',v = are the solutions of the system of equations 
(^-(|l^) for any m at all temperatures. Now let us investigate the behaviour of 
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IRSB system near the possible point of instability. To do this we expand the free 
energy (^) in series near the arbitrary point tc- 

F = Fo + \xl (F,, + rF,,t) + \pl {F^p + rF^pt) + 
+ xipi {F^p + TF^pt) + (F™ + TF^^t) + 
+ xiv {F^x + rF^^t) + Piv {F^p + rF^pt) + 

~l~ '^•^iFxxx ~l~ "^Pi-^ppp ~l~ Fyyy-\- (-^3) 

~l~ 2'^l iPl^xxp ~\~ "vFxxv) H l~2'^l {piFxxp ~\~ "vFxxv) ~\~ 

~l~ 2^1 (■^I'^^PS; ~l~ "^Fppjj^ H i,-^lFvux ~l~ PlFyyp) ~\~ 

+ XiPiVFxpv, 

where Fq - free energy of the RS solution, r = t — tc, and Xi, pi, v are deviations 
from RS solutions {p = p' + pi,x = x' + Xi). Let us look for small solutions of the 
equations for the extremum condition of the functional (0), following the general 
bifurcation theory (see, for example [^). Introducing new variable A = F/{Nkt) 



and using the expressions for Aab from Appendix 1, we write the extremum condition 
in variables x, p, v in the form: 

xiAxx + PiApx - v{m - l)Apx = Bx, (14) 

xiApx + piApp - v{m - l)App = Bp, (15) 

XiApx + PiApp + v[App - 2mD] = B^, (16) 
where D is defined by the expression 

Ayy = -(m - l)[App - 2mD], 



Bx, Bp, By are bilinear in xi, pi, v, r. The equations (14 - |16]) have always the 
trivial solution that corresponds to RS solution. And at the singular point which 
is determined from the condition that the determinant of the linear homogeneous 
system is equal to zero: 



/A A A 

I -'^xx -'^xv 



det 



Apx App Apy I — (1"^) 

\ A^x Ayp Ay 



'■vv 



the solution ceases to be unique. The condition (|T^ can be rewritten in the form 

m (m - 1) {App - 2D) [AxxApp - AjJ = 0. 

Here the coefficients are for RS solution. Using |^, it is easy to see that 

D' — A A — 4^ 

^ -'^xx -'^pp -"-px 
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nowhere turns to zero (see also |^), and the condition 



App -2D = (18) 

defines the temperature of the possible bifurcation: 

1 

t = tc = . 

1,367 

At this point RS solution is 

x'(te) = -0,581; p'(g = 1,449. (19) 

If the condition (pIS]) is fulfilled, then the fourth equation of the system - the 
extremum condition in m - has the form: 

xiAp^ + piApp - v{m - l)App = S„^, (20) 

so in fact it plays a role of the compatibily condition. 



4. To find small solutions of the system (|14| - |Tg, ^y) in the vicinity of tc, we 
have to (see |3^) solve the system of two equations with the determinant not equal 
to zero D' ^ for variables xi and pi, and to define v and m from the bifurcation 
equations. From the system 

XiA^^ + piAp^ = v{m - l)Ap^ + B^, (21) 



XiAp., + piApp = v{m - l)App + Bp (22) 

we obtain 

xi = ^ ^"^^ " ^^ (23) 
n A — n A 

Pi = v{m-1) + , (24) 

so in the main order 

Xi = 0; Pi = v{m — 1). (25) 

Substituting the obtained solutions into the initial equations (|^ - |16|, |2^) (in 
the l.h.s. we need to substitute the full solutions ( [23| ) and (|2^), while for r.h.s. 
( pSD is enough), we obtain the bifurcation equations, that in our case lead simply to 

Bp = By] Bm = By (26) 

or 

TG^ = -v^iG2 + mG3); (27) 
16 

2m — 1 

T Gi = -V—1-G2 + 2m{G2 - G3) + Sm^G^]. (28) 

m 48 



The expressions for the coefficients Gi and their numerical values for RS solution 
are given in the Appendix 1. Solving (pTj) - (|28|), we finally obtain: 



^ = ^^^P^ = 0,427, (29) 

v = vT] = = 3, 13 (30) 

So, we managed to find the small IRSB solution of the system ([^-[T0|) near T^. 
The IRSB solution x = x',p = p' + vrim — 1), f = -ur at m = 0,427 bifurcates 
continuously from the RS solution. 

5. Now let us investigate the stability of IRSB against the further replica symme- 
try breaking (2RSB). Now we divide each former group of replicas with m elements 
into m/irii groups with mi elements. Parameter q^p is assigned pa, if replicas a and 
P belong to the same smallest group (the number of ps is )^ p2 - ii a and [3 

belong to the same "middle" group, but to the different smallest groups (they are 
n(m-mi) |-]-^gj,g^^ g^^^ pq - i{ a and (3 belong to the different groups of m replicas 
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(their number is "■(" "^) ), As usual (see, for example, P3), we deal with fluctua- 



tions inside one group of m replicas and intergroup fluctuations, connected with the 
corresponding overlap of replicas qap- The term in the free energy, corresponding 
to the intergroup replicon mode, which defines the stability of IRSB solution in our 
case, has the form: 



AFi = ]^vl{m - l)(m - mi) (mi - l)(Po - 2Q4 + R^), (31) 



and 



Po = 2-2t^[4:-4:X + p + v- {p + v)% 
= -2t^[2{p + v)-{p + vf - h], 
R, = -2t^[r,-{p + v)% 
where vi wi (see Appendix 2) are defined by the expressions 

Ps = wi - (m - mi)vi, 

P2 = wi + (mi - l)fi. 

The eigenvalue is 

A' = Po - 2Q4 + ^5 = 2 - 2f[4: -Ax- 3{p + v) + 2t^ + rj, 

where 



^3 = / dy'^- , , > 32 
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r^ = Jdy j^^a^^ . (33) 

The other terms in the free energy are given in the Appendix 2 (|38|). With the 
help of RS and IRSB solutions we can obtain the value of A' near tc- Taking into 
account that A'{tc) = and expanding and in series in variables p, v and t, 
using the expressions 

-A = st^aW{W + a^)), (34) 
op 

2t^a'W{3W + 2a')), (36) 



dp 

^''-^""'-Amt^a^W), (37) 



dv dp 

we obtain for the solution (^), ([2^), 

A.' = QtlmTv{l - 2tlGA). 
Here we take 6*4 for RS solution: 

G4 = {W{a^ + 2a=^ - _ 2a)) = 0, 27, 

so finally 

A' = 2, 97r, 

that means the stability of IRSB solution in this vicinity at T < Tc and mi < m < 1. 

6. In conclusion, we have investigated the behaviour of the model (|l]) near the 
instability point tc of its RS solution. We have shown that in the vicinity of tc there 
is a more favorable solution which corresponds to the one-step replica symmetry 
breaking. We have also shown, that the new solution is stable relative to further 
replica symmetry breaking, i.e. the full Parisi scheme does not work in this case, 
at least in its classical variant. It is worth to notice that in spite of absence of 
reflection symmetry, IRSB solution bifurcates continuously from the RS solution, 
that contradicts the idea of division in two classes. 

However, one can await that at further decrease of the temperature, our IRSB 
solution will become unstable. For complete investigation of this case it is necessary 
to solve the system of equations ([THIDD for all temperatures and use the obtained 
solution for the definition of the sign of Ai everywhere. May be, the low-temperature 
form of the IRSB solution will explain the behaviour of experimental data at 



ultralow temperatures T ^ 0, 05 where the deviations from the RS solution are 
noticed, i.e. the linear dependence of the specific heat on T changes to the quadratic 
one. 
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1 Appendix 1 

Let us take into account that f = for the RS solution. So \E' does not depend on 
and the internal integrals in (H - |^) can be taken directly. Let us denote: 



a, 



In these notations 



t=tc, x=xo,p=po, y=0 
W = a^ + a-2, 
(...) = / dy^... . 

Ap^ = --{Wa), 

App = -l + t^{W{W + 2a^)), 
Avx = —{m — l)Apx, 
Apv = —{m — l)App, 
A^y = -(m - l)[App - 2mD], 
D = t^{Wa^), 
App-2D = -l + t^W^), 
At the bifurcation point t^. = x' = —0, 581, p' = 1, 449 and we have for the 



RS solution 



4,132, 



(Wa) = 0,6035, 
(Wa^) = -0,997, 

(Wa^) = 1,232, 
{Wa^) = -1,981. 



The coefficients Br^, Bp, By contain the third and B^ - the fourth derivatives of 
A at the point and their numerical values are obtained from the RS solution. 
With their help we get: 
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Gi = -t + t^W[{x'{tc) - 2)(2a3 + - 3a - 2) + 
+p'{Qi-10a^ + 18a^ + 15a^ + 19a - 6)]) = -0,4, 

G2 = 4(l^[4a^ + 8a^ - 3a^ -7a- 2]) = 11, 8, 
G3 = 4(W^[-10a^ - 2Qa^ + 12a^ + 22a + 10]) = 4. 



2 Appendix 2 

4AF2 = p^m[-2Pi + 4mQi - 4(m - 1)Q3 - Qm^Ri + 8m(m - l)i?3 - 2(m - l)2i?7] + 
+ w^im - l)[2Po + 4Q4(m - 2) + i?5(m - 2)(m - 3) - i?6m(m - 1)]+ 
+ PoWiAm{m - 1)[-2Q2 + mi?2 - (m - 2)i?4]. 



Here the coefficients P, Q, R are for IRSB solution and are as foUows: 



Pi = 2-2i^[4-4x+p-p^], 

Qi = -2i'[2p-p2-tni], 

Q2 ^ -2t^[2p - p{p + v) - hi], 
= -2t^[2{p + v) - p^ - t2i], 
Ri = -2t^[rnii - P% 
R^^-2f[r2u-pip + v)], 

R3 ^ -2f[r2ii - p% 
R4^-2f[r3i-p{p + v)i 
R6^-2t''[r22-{p + v)% 
R'r^-2e[r22-p\ 



(38) 



where 




(41) 



(43) 



(39) 



(40) 



(42) 
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